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ABSTRACT 

Avalanche photodiodes (APD) are studied for use as photon counting detec- 
tors in space-borne lidars. Non break-down APD photon counters, in which the 
APDs are biased below the break-down point, are shown to outperform conven- 
tional APD photon counters biased above the break-down point and to outper- 
form APDs in analog mode when the received optical signal is extremely weak. 
Nonbreak-down APD photon counters were shown experimentally to achieve an 
effective photon counting quantum efficiency of 5.0% at X = 820 nm with a 
dead time of 15 ns and a dark count rate of 7000/s which agreed with the 
theoretically predicted values. The interarrival times of the counts followed an 
exponential distribution and the counting statistics appeared to follow a Poisson 
distribution with no after pulsing. It is predicted that the effective photon 
counting quantum efficiency can be improved to 18.7% at X =820 nm and 
1.46% at X = 1060 nm with a dead time of a few nanoseconds by using more 
advanced commercially available electronic components. 
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1. Introduction 

Space-borne lidar systems have become an important tool for earth science 
and space explorations. Diode laser pumped NdfYAG lasers are among the most 
attractive laser sources for space-borne lidars because of their small size, high 
output power, high efficiency, high optical beam quality, long lifetime, and 
inherent system redundancy. Silicon avalanche photodiodes (APD) can provide 
the highest receiver sensitivity among other types of photodetectors at the fun- 
damental wavelength of Nd:YAG lasers (X=1060nm). APDs may be used in 
either analog or photon counting (Geiger) modes. We report here a detailed 
study of the sensitivities of lidar receivers which use APD photodetectors, espe- 
cially APDs operated in the photon counting mode, when biased slightly below 
the break-down voltage. APDs under this operation condition are shown to 
always outperform APD photon counting receivers biased above the break-down 
voltage and to outperform APDs in analog mode when the incident light level is 

extremely low. 

Lidars in space have mainly three types of scientific applications: cloud and 
aerosol ranging [l], atmospheric back light scattering measurement [2], and 
spacecraft to ground altitude ranging [3]. All these applications require the lidar 
receivers to measure the reflected light energy integrated over a short time inter- 
val. The lidar receivers have to be extremely sensitive and have wide dynamic 
ranges because of the great distances from the host spacecraft to targets and the 
possible variation in the reflectivities of targets. Silicon APDs are the photo- 
detectors of the choice at X=1060nm because of their relatively high quantum 
efficiencies which are typically 6% [4], but can be as high as 40% [3, p. 468]. As 
a comparison, the quantum efficiencies of typical photomultiplier tubes (PMT) 
are only about 0.1% at the this wavelength [5]. 
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APD’s are normally operated in either analog or avalanche break-down pho- 
ton counting modes. In analog mode, an APD is biased slightly below the 
break-down voltage such that the average output current is linear in the incident 
optical power and equal to the primary photocurrent multiplied by an internal 
avalanche gain. In the avalanche break-down photon counting mode, an APD is 
biased above the break-down voltage such that individual detected photons 
trigger an avalanche break-down which results in a strong output current spike 
that can then be individually counted. The performance of APDs in analog 
mode are limited by the APD excess noise due to the randomness of the APD 
gain mechanism. On the other hand, APDs in avalanche break-down photon 
counting mode are subject to after pulsing and a relatively long dead time due 
to the slow recovery process after each avalanche break-down. 

APDs may also be operated in photon counting mode with the bias voltage 
below the break-down point. The average APD gain in such a photon counter 
should be set to its highest stable value, which may require APD temperature 
control and a highly stable bias voltage supply. Each photon absorption will 
result in a photocurrent pulse but the amplitude of these photocurrent pulses are 
much smaller than those output by APDs biased above the break-down point. 
Consequently, a high gain low noise amplifier and discriminator are required to 
detect these relatively weak photocurrent pulses. The major advantages of this 
type of photon counter over the break-down APD photon counters are the elimi- 
nation of both after pulsing and the slow recovery process following each 
avalanche break-down. As a result, the counter dead times is limited only by 
the electrical bandwidth of the APD and the subsequent circuits. The major 
obstacle toward implementation of this type of APD photon counter has been 
the lack of wide band and low noise APD preamplifiers. Kikuchi and et. al. 
(1985) first implemented and tested such a nonbreak-down APD photon counter, 
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in which the APD and preamplifier were held at 77 ° IC in order to reduce ther- 
mal noise [6]. Lightstone and McIntyre (1988) gave a partial characterization of 
this type of APD photon counter along with some measurement data at room 
temperature [7]. However, the dead time of their experimental photon counter 
was rather long (>0.5 /is) because the electrical bandwidth had to be kept rela- 
tively small in order to find a preamplifier with sufficiently low thermal noise. 

We have developed a complete theoretical model to characterize this type of 
nonbreak-down APD photon counter. An experimental system was built and 
tested. The measured data agreed well with the theoretically predicted values. 
The measured photon counts appeared to have a Poisson distribution and the 
interarrival times fit an exponential distribution, as expected. The dead time 
was measured to be about 15 ns, which is the shortest yet reported for APD 
photon counters operated in either nonbreak-down or break-down mode. The 
overall photon detection efficiency was about o.0% at X=820nm, which projects 
to about 0.39% at X=1060nm, based on the typical spectral response curve of 
silicon APDs [4]. It is predicted that the overall photon detection efficiency can 
be improved to as high as 18.7% at X=820nm and 1.46% at X=1060nm with a 
dead time of a few nanoseconds by using more advanced electronic components. 

The rest of this report is organized as follows. Sections 2 and 3 give brief 
reviews of APDs in analog and break-down photon counting modes, respectively. 
Section 4 gives a detailed description and performance analysis of APD photon 
counters under nonbreak-down operation. The experiments and measurement 
data are presented in Section 5. Section 6 contains conclusions along with some 
predicted performance of APDs in nonbreak-down photon counting mode when 
usin 0- more advanced electro-optical and electronics devices. 
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2. APDs in Analog Mode Operation 

When an APD is biased below the break-down voltage, the average output 
photocurrent is proportional to the incident optical power. The photocurrent in 
a lidar receiver is usually first integrated and then fed into a microprocessor 
through an A-D converter. The number of photons received over the integration 
time can be estimated by dividing the integrated photoelectron charge by the 
electron charge and the average APD gain. The advantage of APDs over ordi- 
nary photodiodes is that the primary photoelectron current in an APD is multi- 
plied by an internal multiplication gain, as in a photomultiplier tube. Since in 
practice there is always circuit thermal noise present which appears as additive 
noise in the total APD output photocurrent, the use of an APD can increase the 
level of the optically generated signal and consequently improve the receiver sen- 
sitivities. However, use of an APD also causes an excess noise in the APD out- 
put due to the randomness of the APD gain. The excess noise is usually charac- 
terized by an excess noise factor defined as F =E {m 2 }/(E {m}) 2 with m the ran- 
dom multiplication gain of the APD. The excess noise factor is related to the 

APD device parameters by [8] 

F =k eff G +(2~)(1 -K ff ) (1) 

where £ e yy is the ratio of the ionization coefficients of holes and electrons in the 
multiplication region of the APD, G is the average APD gain which is controlled 
by the bias voltage. The smallest value of the excess noise factor is two for ideal 
APDs with k e f f — *0 and G — kx>. 

In practice, the output signal from an APD is still too small to be con- 
veniently used and it must be further amplified for subsequent signal processing. 
Each amplifier adds thermal noise to the output but the dominate thermal noise 
source is from the first stage amplifier which is often called the preamplifier. 
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There are in general two classes of preamplifiers: high input impedance amplifiers 
and transimpedance amplifiers, as shown in Figures 1 and 2. The former is a 
high input impedance linear voltage amplifier preceded by an APD load resistor. 
The latter is a feedback amplifier with the feedback resistor as the APD load 
resistor. 

Figure 1 also shows an equivalent circuit diagram of a high input 
impedance preamplifier. The bandwidth of the preamplifier can be approxi- 
mated as inversely proportional to the time constant of the APD load resistance 
and the total shunt capacitance at the front end of the amplifier. The input 
shunt capacitance is very difficult to reduce below 1 pF unless the APD and the 
amplifier are integrated into one chip. The APD load resistance should be as 
large as possible because the thermal noise current generated by the resistor is 
inversely proportional to the feedback resistance. In order to achieve low noise 
and wide bandwidth operation, an equalizer is usually used after the preamplifier 
to enhance the high frequence response of the entire circuit. The equalizer 
should, in fact, cancel the pole formed by the input capacitance and the APD 
load resistance. Since the input capacitance is parasitic and drifts with tempera- 
ture, a good equalizer is hard to build and the entire circuit may even become 
unstable. The equalizer also enhances the high frequency noise of the amplifier 
which becomes dominant as the receiver bandwidth exceeds a certain value. 
Furthermore, the APD load resistor at the input of the amplifier loads down the 
bias voltage across the photodetector which may become a major factor in limit- 
ing the dynamic range of the receiver. Therefore, high input impedance 
preamplifiers are not used very often in practice, though they can potentially 
achieve a higher receiver sensitivity than the other type of preamplifiers [9], 

A transimpedance preamplifier is easier to build and more stable in opera- 
tion. An equivalent circuit of a transimpedance preamplifier is also shown in 
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Figure 2. The transfer function of the preamplifier can be derived as 


Zt (w) — 


1 (cj) 

A[oS) Z { (oj) 


aZ l^- T+J^C i’ " ( 2 ) 

where Z f {oJ) and Z,(w) are the total input and feedback impedances, A (w) is 
the open loop amplifier gain, and R y and C y are the feedback resistance and its 
parasitic shunt capacitance. As shown in (2), the 3dB bandwidth of a transim- 
pedance amplifier can be approximated as (27T RjCj') * Hz. Since the parasitic 
capacitance across the feedback resistance can be made very small with modern 
integrated circuit technology, the 3dB bandwidth of a transimpedance amplifier 
is usually much larger than that of a high input impedance preamplifier with the 
same value of the APD load resistance. In fact, there is usually no need for an 
equalizer following the preamplifier. Furthermore, a transimpedance amplifier 
has little loading effect on the photodetector bias voltage and has a wider 
dynamic range than a high input impedance preamplifier. Therefore a transim- 
pedance preamplifier is preferred from the overall system performance point of 
view. The performance analysis in the subsequent sections always assumes that 
transimpedance preamplifiers are used. 

It is convenient to define the received signal as the number of photoelectrons 
accumulated over the integration time of the hdar receiver. The mean and vari- 
ance of the signal can be computed as [10] 


X - a?g-(P,+ P it )T + 

hf 


GI b r | I s t 
Q ? 


( 3 ) 


and 
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o? = FG 2 ^-(P, + P„ )t + FG 2 ^- + -y- + <4n P (4) 

where t]q is the quantum efficiency of the APD, hf is the photon energy, P 3 and 
P bg are the received signal and background optical power, r is the integration 
time, I a and I b are the .APD surface and bulk leakage currents, q is the electron 
charge, and is the variance of the total amplifier noise equivalent at the 

input of the preamplifier integrated over t. Usually, high speed transimpedance 
preamplifiers have an FET front end and the variance of the total amplifier 
noise can be expressed as [1 1] [12] 


"f- + I J + 

If 9 q“Rr 2 7 


where K is Boltzmann’s constant, T a is the temperature in Kelvin, R f is the 
feedback resistance, I g is the gate leakage current of the FET, E { is the one 
sided power spectrum of the series input voltage noise of the amplifier and is 
given by Ej == 4KT a T'/ q m with q m the transconductance of the FET and T a 
numerical factor close to unity, B c is the cutoff frequency of the integrator, and 
C { is the shunt capacitance at the input of amplifier. Theoretically, the fre- 
quency response of an ideal integrator is a sine function, and therefore, the 
bandwidth required to implement it must be infinite. In practice, the bandwidth 
is limited to several times the main lobe of the sine function. It has been shown 
experimentally [13] that a cutoff frequency of B c >2(1/t) is sufficient for the 
degradation in performance to be negligible. 

The first two terms in (5) are proportional to r and can be treated as hav- 
ing a constant spectral noise current density (white noise). The last term in (5) 
corresponds to a spectral noise current density which is proportional to the 
square of the frequency. Therefore, it is important to keep the parasitic shunt 
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capacitance at the input of the amplifier as low as possible. For most commer- 
cial preamplifiers which have moderate feedback resistance (~Kf2), the power 
spectrum of the total amplifier noise can be assumed to be constant (white) at 
least up to a few hundred megahertz. Under this condition, an equivalent noise 
temperature, T e ( ° K), can be defined such that the total amplifier noise is equal 
to the thermal noise generated by the feedback resistance at temperature T e . 
The total amplifier noise can then be written, to a good approximation, as 


amp 


2 KT e r 

7r } 


( 6 ) 


If we define the useful signal as x s =G [ t)q / hf )P s t, the signal to noise ratio 
(SNR) of the APD output in analog mode can be written as 


-2 


SNR an = 




FG 2 


0q 

hf 


J b T I s t 2KT e r 

<P a + P kg )r + FG 2 ±- + — + 


(7) 


q 2 Rj 


There exists an optimal value for the average APD gain which maximizes the 
SNR an . At the optimal APD gain, the excess noise should be much larger than 
the amplifier noise. For comparison purposes, we can approximate the SNR an as 

Sm„ -» jVqSNR ql (S) 

where SNR QL is the quantum limited SNR assuming unity quantum efficiency 
and no amplifier and background radiation noises. Therefore, roughly speaking, 
the SNR an is a factor of F smaller than that achieved with a photodiode with 
the same quantum efficiency but no thermal noise and no excess noise. In other 
words, use of an APD makes the amplifier noise negligible but also causes the 
signal to noise ratio to decrease roughly by a factor of F. Nevertheless, APDs 
can significantly improve receiver SNR at X=1064nm as compared to PMTs 
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which have little excess noise (F~l) but a much lower quantum efficiency. For 
example, if we assume F=8.5, 30%, and ? ?Q PMr = 0-l%, the use of an APD 

can achieve a receiver SNR which is 35 times that of a PMT, except when the 
optical signal is so weak that the required optimal APD gain is too high to be 
obtained in practice. 


3. APDs under Break-Down Photon Counting Operation 

An APD can be biased several volts above its break-down point to achieve 
photon counting operation [4]. Avalanche break-downs are triggered by both 
absorbed photons and thermally generated hole-electron pairs inside the device. 
Each avalanche break-down results in a large current spike which can be easily 
detected by a discriminator. Figure 3 shows a diagram of a typical APD break- 
down photon counting circuit [4]. 

Since not all the received photons can trigger avalanche breakdowns, the 
average number of photons counted by the discriminator, n c , is always less than 
the average number of absorbed photons, n. The effective photon counting 
quantum efficiency can be written as 

VQ, f} =VQ^=- = V Q Ptri g ( 9 ) 

ft 

where r] Q is the quantum efficiency of the APD, P trig is the probability that an 
absorbed photon triggers an avalanche break-down which results in a discrimina- 
tor threshold crossing. The average number of absorbed photons, n, over a 
counting interval, T, is related to the received optical power, P 0 , by 

Z = r, Q P 0 T/hf. (10) 

The number of photons being counted, n, should be a Poisson random variable 
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with mean and variance equal to n c . The signal to noise ratio is given by 

SNR cl - ~n c = P, ri , n (11) 

0 var{n } * 

Considering the background noise counts and dark counts, the signal to noise 
ratio should be 


SNR ct = 


{£t rigjhf_ 

P trig{ n a^ r n bg ) "F n dark 


( 12 ) 


where n a and n bg are the average number of absorbed signal and background 
photons, and n dark is the average number of total dark counts caused by the 
APD dark current and circuit thermal noise. 

The probability of triggering, P trig , increases with the bias voltage. How- 
ever, the number of dark counts triggered by thermally generated hole-electron 
pairs also increases with the bias voltage. The value of Ptrig is typically 7% 
with 15000/sec dark counts at 22 °C when the APD is biased approximately two 
volts above the break-down point [4]. Although higher values of P trig (50-80%) 
can be achieved at even higher bias voltages [14], they cannot be used in lidar 
receivers because of the excessively high dark count rates and other nonlinear 
effects. The dark counts can be reduced by cooling the APD. However, the 
quantum efficiency of the APD may also decrease with temperature over certain 
wavelength ranges. For example, for a RCA30902S silicon APD, the dark count 
rate decreases from 15000/sec at 20 0 C to 300/sec at -25 ° C, but the quantum 
efficiency also decreases by about a factor of two at \=1064nm [4]. 

A major disadvantage of break-down type APD photon counters is the 
existence of so called after pulsing due to excess avalanche break-downs which 
appear randomly after each avalanche break-down. The cause of after pulsing is 
attributed to the presence of trapped charges at dislocations and impurities in 
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the high voltage region of the APD after each avalanche break-down [7]. After 
pulsings can occur randomly at nanoseconds to kiloseconds after each avalanche 
break-down but they are mainly distributed within the first a few milliseconds 
[4] [14]. The probability of after pulsing within 1 microsecond to 60 seconds is 
typically 2-15% [4] and increases with the APD bias voltage. After pulsing is a 
nonlinear effect and has not been well understood except for some simplified 
mathematical modeling [15]. In lidar applications, after pulses appear as ran- 
dom echos from the target which reduces the SNR ct or causes confusion in the 
measurements. Multiple laser firings at one target spot may be required to aver- 
age out the random echos. In photon correlation spectroscopy measurements, the 
decay times of the after pulses form a systematic measurement error which some- 
times can be overwhelming. 

The dead time of a break-down APD photon counter is mainly determined 
by the recovery process of the APD after each avalanche break-down. Every 
time a break-down occurs, the APD has to be quenched to prevent any damage 
to the device. There are two types of quenching circuits: passive quenching and 
active quenching. One example of passive quenching is shown in Figure 3. 
When a break-down occurs, the voltage drop across the resistors in series with 
the APD effectively reduces the APD bias voltage to below the break-down 
point. The APD bias voltage is then gradually restored as the APD shunt capa- 
citance gets recharged. Since the resistors have to be relatively large to limit the 
maximum APD current, the recovery process is relatively slow. The typical dead 
time of this type of photon counter is about one microsecond [4] [16]. An active 
quenching circuit can significantly shorten the dead time by momentarily reduc- 
ing the resistance in series with the APD upon each break-down so that the APD 
shunt capacitance can be charged with a much shorter time constant than that 
of a passive quenching circuit. Brown and et. al. [17] have reported an active 
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quenching circuit which achieved a dead time of a few tens of nanoseconds. 

Because of the existence of a finite dead time, photons which arrive within 
the dead time cannot be resolved and the number of detected photons is no 
longer a Poisson random variable. Cantor and Teich [18] have shown that the 
probability that n photons are counted over a T second interval in response to 
n average absorbed photons for a given dead time, T d , can be written as 


Pd( n ' \ n > T d ) = 


E 

Jt-0 


[n(l ~n'f 
~k\ 




*-0 


WH*'-l)y}' - 
e 


»[!-(«'■ -ll-yl 


n ' < • (13) 

T d 


The maximum likelihood (ML) estimator, n ML , of the average number of actu- 
ally absorbed photons, n, given n detected photons under this condition is 
given by [19] 


HfL 


n 


1 “ n'i T d / T ) 


(14) 


As Tj =0, Eq. (13) reduces to a Poisson distribution and Eq. (14) reduces 
t° fi ML = n = n which is the actual number of photons being counted. The 
mean and variance of the ML estimator are given by 

E {n} = va,r {n} = n . (13) 

The signal to noise ratio under this condition is equal to 

SNR 0 - = n. (16) 

var (n } 

When T d ± 0, the ML estimator given by (14) become biased and has a larger 
variance, i.e., (E{n}-n )± 0 and As a result, the signal 
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to noise ratio becomes smaller, compared to the ideal case, i.e. 

SNR (») = - = SNR 0 . 

var {n } var {n } 


( 17 ) 


Figures 4 and 5 show the normalized bias and signal to noise ratio as a function 
of average number of absorbed photons. It is shown that {r d /T) < 0.01 is 
required in order to obtain a reasonably good measurement (i.e. bias~0 and 
SM?(n)~n. The dead times of avalanche break-down APD photon counters are 
usually too long for lidar receivers. 


4. APDs under Nonbreak-down Photon Counting Operation 

APDs may be used to count single photons while being biased below the 
break-down point as in analog mode [6] [7]. The photocurrent pulses resulting 
from a single photon absorption can be detected via a discriminator threshold 
crossing with a relatively high probability as long as the APD guiu is sufficiently 
large and the circuit thermal noise is relatively low. There is no after pulsing 
since the APD is operating in its linear regime. The dead time can also be much 
shorter than that of a break-down type APD photon counter since the APD bias 
voltage does not need to be restored after each photon detection. The receiver 
can be operated in both photon counting and analog modes simultaneously by 
diverting part of the APD preamplifier output to an analog channel. The analog 
channel can operate linearly after the photon counting channel is saturated due 
to a strong input light intensity. Therefore, the dynamic range of the entire 
receiver can be greatly extended. 


- 15 - 


4.1. Principle of Operation 

The circuit used for this type of photon counter is shown in Figure 6. The 
APD and preamplifier are exactly the same as those in the analog mode except 
that the average APD gain is set to the highest achievable value but not neces- 
sarily the value which optimizes the output SNR an . This requires some form of 
APD temperature control or compensation and a very stable bias voltage supply. 
The output of the preamplifier is compared against a threshold with the use of a 
comparator or a discriminator whose output is connected to an electrical pulse 
counter. The threshold level of the discriminator has to be carefully set in order 
to achieve a high photon detection efficiency while maintaining a tolerable noise 
count rate. As shown in Figure 6, the dead time of this type of photon counter 
is limited by the photocurrent pulse width which is about equal to the reciprocal 
of the bandwidth of the APD and preamplifier. In practice, the dead time is 
longer than the pulse width because of the limited time resolution of the discrim- 
inator and the electrical pulse counter. Since the APD gain is random, the pulse 
width is a random variable and so is the dead time. The actual dead time distri- 
bution usually has to be determined experimentally. 

The number of detected signal photon counts and the number of dark 
counts are both Poisson random variables. The signal to noise ratio can be com- 
puted using Eq. (12) by substituting the threshold crossing (discrimination) pro- 
bability, P di3C , for the triggering probability, P trig . 

4.2. Effective Quantum Efficiency 

Since the APD gain is random, not all the photon absorptions result in pho- 
tocurrent pulses whose peak amplitudes are large enough to cross the discrimina- 
tor threshold. The effective quantum efficiency can be written as 
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_ E{photons_ counted^ „ 

^Q'/f E {photons absorbed } ^ d,sc 


(18) 


where rjq again is the quantum efficiency of the APD and P disc is the probabil- 
ity that a photocurrent pulse in response to a photon absorption is counted. 
The computation of P disc is given as follows. 

Ideally, the photocurrent pulse output from an APD in response to a pho- 
ton absorption should be an impulse and the preamplifier under this condition 
can be simplified equivalently to consist of a resistor and a capacitor in parallel. 
The pulse shape output from the preamplifier is an exponential function as 
shown in Figure 7 (dotted line). The average area under the pulse shape function 
is mq with m the instantaneous APD gain and q the electron charge. The peak 
amplitude of the pulse is given by i phtnax =mq /RC=2nGqB 3dB with B 3dB (Hz) 
the 3 dB bandwidth of the preamplifier. In practice, the APD output cannot be 
true impulses and the preamplifier is always bandwidth limited. As a result, the 
shapes of the actual photocurrent pulses output from the APD are not true 
exponential functions but have a finite rise time, a round peak, and a slower fal- 
ling edge, as shown in Figure 7 (solid line). Empirically, we choose 


ipfimaz ( 1 ^) 

which seems to yield results consistent with the measurements. Equation (19) 
may have to be modified if the frequency response of the preamplifier or the sub- 
sequent circuit is significantly different from that of a RC lowpass filter. 

A photon is detected whenever the peak amplitude of a photocurrent pulse 
exceeds the threshold current, I thre . The probability of discrimination, P d ; sc , is 
then given by 


Pdisc = Prob (i phmaz > Ithre ) = Prob ( m ^ M thre ) 


CO 

- £ p ( m ) 

m — Mikre 


( 20 ) 
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where M thre is the equivalent threshold crossing APD gain given by 


Mt lire ~ 


^thre 

^Q^3dB 


( 21 ) 


and P{m) is the the probability of the APD gain. McIntyre and Conradi 
[20] [21] have shown that the probability of the APD gain can be written as 
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with r( ) the GAMMA function. The equivalent quantum efficiency of this type 
of photon counter can be evaluated by substituting (21) and (22) into the right 
hand side of (20). 

Figure 8 shows qualitatively the distribution of the APD gain for two 
different values of the ionization ratios k e jj- [20]. The probability of discrimina- 
tion P disc is equal to the shaded area under the curves. It is seen that the 
smaller the k e j f , the larger the P digc . Figure 9 plots P disc as a function of the 
normalized equivalent threshold crossing APD gain, M^/G, for G=1200 and 


k e f j =0.010. 


4.3. Dark Counts 

There are two sources of dark counts: the APD dark current and the circuit 
thermal noise. The APD dark current consists of surface leakage current and 
bulk leakage current. The former does not pass through the high gain region of 
the APD and hardly causes any dark counts. The latter does pass through the 
high gain region of the APD and acts as a source of background radiation. The 
number of dark counts due to a bulk leakage current, , can be modeled as a 
Poisson random variable with the mean equal to 
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( 23 ) 

where P disc is the average probability that an electron of the APD bulk leakage 
current causes a discriminator threshold crossing. The value of P disc should be 
smaller than that for a primary photo electron, i.e. P disc < P d { sc , since part of 
the APD bulk leakage current is generated inside the APD high field region and 
is not multiplied by the full APD gain. We have found experimentally that 
P ddac ~ °- 17P disc f° r a RCA30902S silicon APD which had a bulk leakage current 
of / 6 ~ 0.1 pA 

The number of dark counts due to circuit thermal noise should also follow a 
Poisson distribution if the correlation time of the noise is less than the dead time 
of the photon counter. The average number of counts due to circuit noise, 
n dark2 , o y er a counting interval of T seconds can be written as [22] 

n dM - Pf . ( 24 > 


where At is the average threshold crossing time or pulse width which is about 
equal to the reciprocal of the band width, and P f a is the threshold crossing pro- 
bability due to the circuit thermal noise at a fixed sampling time. The circuit 
noise current can be modeled as a zero mean Gaussian random process and the 
variance can be written similarly to (6) as 


~2 = 
w lamp 


4KT, 
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(25) 


where B n is the noise bandwidth. If the transimpedance preamplifier can be 
assumed equivalently as an ideal RC circuit, the noise bandwidth is 1.57 times 
the 3dB bandwidth [23]. The threshold crossing probability can be expressed as 

oo 

p. = f 1 - e 2 ° lm ’ dx . 

V 27rcr /amp 
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There is always a trade-off between the equivalent quantum efficiency and 
circuit noise counts. Both P d { sc an d n dark2 decrease as the threshold I tkre 
increases. Figure 10 shows qualitatively the average number of detected photons, 
nP di S ci and noise counts, n dar k2 , as a function of the threshold current /^ re . For 
best performance, the threshold level should be set at the point where the curve 
of n dark2 just starts to level off and the total number of dark counts is much less 
than that of the photon counts, i.e. n darkl + n dark2 « n -P disc . If the value of 
the number of absorbed photons, n , is completely unknown and needs to be 
determined through the measurement, the threshold level should be set such that 
the dark counts due to circuit noise are much less than those due to the APD 
bulk leakage current, i.e. n dark 2 « n darki • Usually, the threshold is set to be a 
few times the standard deviation of the circuit noise, i.e. I thre = / Xc Iamp with 
/ a multiplication factor. The equivalent threshold crossing gain in (21) can be 
expressed as 


AKT, 


f )] 1/2 


^thre 


R 


MPzdB 


(27) 


The probability of discrimination and the equivalent quantum efficiency can be 
evaluated by substituting (27) and (22) into (20) and (18). If the photon counter 
is operating properly, the numbers of detected signal and background photons, 
the dark counts due to APD leakage current, and the dark counts due to circuit 
thermal noise should all follow Poisson distributions, and so should the total 
number of counts. 

As seen in (27), the equivalent APD threshold crossing gain M thre is 
inversely proportional to the square root of the feedback resistance Rf and the 
3dB bandwidth of the preamplifier. Since the probability of discrimination P disc 
increases as M thre decreases, it is important to keep the product of the feedback 


- 20 - 


resistance and the 3dB bandwidth of the transimpedance preamplifier as large as 
possible in order to achieve a high photon detection probability. 


5. Experiments with Nonbreak-down APD Photon Counters 

An experimental nonbreak-down APD photon counter was built and the 
performance was measured to verify the theoretical model developed in the previ- 
ous section. 

5.1. Experimental Setup 

Figure 11 shows the optical setup of the experiment. The light source con- 
sisted of an ordinary incandescent light bulb, a diffuser, a pinhole, and a col- 
limating lens. The light bulb was powered by a stable DC voltage supply and 
the light emitted can be considered to have an almost perfectly constant inten- 
sity over a wide spectral range. The light beam was first attenuated and then 
filtered through two interference filters which were centered at X=821nm and 
had a combined bandwidth of AX=7nm. The received optical power was meas- 
ured directly by substituting an optical power meter sensor head for the APD 
which could be moved aside. The experiment was conducted in total darkness 
and background radiation could be neglected. A high speed universal counter 
(Stanford Research SR620, maximum count rate 200X10 6 /sec ) was used to 
count the pulses output from the discriminator. The counter also gave the 
mean, the variance, and histograms of the measured counts. 

The details of the electronics are shown in Figure 12. The APD used was a 
RCA 30902S silicon APD with Jfc e// = 0.010, / 6 ~0.1 pA , and I s = 12 nA [13]. 
The APD preamplifier consisted of a GigaBit Logic 16G071 transimpedance 
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amplifier (DC-700 MHz, Rj= 830 Cl [24]). The signal was further amplified by 
another two amplifiers to a level which was appropriate (~200 mV) for the 
discriminator. The attenuators shown in Figure 12 were used to adjust the 
amplifier gains since there were no other amplifiers in our lab which had just the 
right gain and dynamic range. A low pass filter was used to block out the high 
frequency noise and the cutoff frequency (200 MHz) was chosen to correspond to 
the maximum speed of the discriminator and the counter. The discriminator 
consisted of a high speed comparator (Motorola MC1650, MECLIII series one bit 
A/D converter) and the resolution of the comparator was smaller than 20 mV 
according to the data sheet. 

A schematic circuit diagram of the APD/preamplifier DC voltage supply is 
shown in Figure 13. The high voltage APD bias supply was generated by a 
DC-DC converter (Analog Modules, Model 521) which had peak output ripples 
of <0.003% of the average output DC voltage. 

5.2. Circuit Thermal Noise 

The total circuit noise at the input of the comparator was measured with a 
RF power meter (HP-435B with 8482A power sensor head) while biasing the 
APD about 100 V below its break-down point so that the APD gain was essen- 
tially zero and its noise could be neglected. The total RF noise power measured 
was 2.5 fxW over a bandwidth of 0.1-100 MHz. The net gain from the output of 
the preamplifier to the input of the discriminator was measured to be 46.5±0.5 
dB from 0.1 to 200 MHz with the -3 dB point at 220 MHz. The resultant 
equivalent noise current at the front end of the transimpedance amplifier was 6.4 
pA/\ZHz, which was very close to that given by the data sheet (5.7-6. 3 
pA/VHz) [24]. Figure 14 shows the power spectrum of the noise measured by a 


- 22 - 


RF spectrum analyzer (HP-8590A). The spikes in the spectrum were due to 
electro-magnetic interferences (EMI). The average noise density up to 100 MHz 
was about -107 dBm/Hz (2.0 /AV/100MHz) according to Figure 14 which was 
close to that measured by the RF power meter (2.5 /fW/100MHz) and the 
discrepancy was believed to be caused by the spikes in the power spectrum. The 
noise intensity rose by a few dB over the frequency range of 100-200 MHz as 
shown in Figure 14. The total noise power measured with the RF power meter 
at the input of the comparator was about 7.0 fSW (18.7 mV rms) over the entire 
bandwidth. 

The average noise counts due to the circuit thermal noise were first meas- 
ured as a function of the absolute value of the threshold level applied at the 
comparator. The result is shown in Figure 15. The APD was biased well below 
its break-down point so that the contribution from the APD noise could be 
neglected. The counting interval was set to one second in this measurement and 
each point in Figure 15 was based on the average of 100 measurements (sample 
size). The threshold voltage was measured at the negative input terminal of the 
comparator with the use of a digital multimeter. The actual threshold voltage 
was negative since the transimpedance APD preamplifier was an inverting 
amplifier. It is noticed that the effective threshold level defined in the previous 
section was higher than that measured at the comparator since the comparator 
had only limited gain and pulses at the positive input terminal had to exceed the 
threshold by a certain amount in order to trigger the subsequent counter. For 
example, the measured threshold voltage was -47.1 mV at a dark count rate of 
900 cts/sec and the effective threshold level which corresponded to the same dark 
count rate should be about 4.5 times the total measured rms circuit thermal 
noise voltage at the input of the comparator according to (24)-(26), i.e. V e jf th = 
4 > 5X(— 18.7m V) = —84.1 mV. Therefore, the difference between the measured 
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and effective threshold levels was about 37 mV in this case. 

5.3. Measurements of Photon Counts and Dark Counts 

The APD in this measurement was biased as close to the break-down vol- 
tage (~252 volt) as possible without an avalanche break-down occurring at any 
time. An oscilloscope was used to monitor the APD preamplifier output to detect 
any possible avalanche break-downs which appeared as abnormally large spikes 
in the output waveform. The actual average APD gain could not be measured 
directly because the APD was under CW illumination and the amplifiers were 
AC coupled. However, we had measured the average gain of another APD of the 
same model number under pulsed illumination and the highest average gain 
achieved was found to be 1200-1300. Figure 16 shows the average waveforms of 
input pulses (negative going) and output pulses (positive going) of the compara- 
tor. It is noticed that the measured average pulse width (5.857 ns) was only 
approximated since it depended on the triggering level of the oscilloscope. 

The received optical power was adjusted to 1.0 pW. The numbers of total 
counts and dark counts were measured as a function of the absolute value of the 
threshold voltage. The means and the standard deviations (error bar) of the 
measured counts are plotted in Figure 17. The dark counts were measured while 
blocking the incident light beam. The counting interval was set to 1.0 ms and 
the sample size was 10,000. The means and standard deviations of the thermal 
noise counts for the same counting interval and sample size are also shown in 
Figure 17 (dot-dashed curve). The number of photon counts increased relatively 
slowly with the threshold level as compared to the dark counts which included 
the circuit thermal noise counts. At low threshold levels, the circuit thermal 
noise counts became overwhelming and the standard deviations became 
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abnormally large, as shown in Figure 17. The optimal threshold should be at 
the "knee" of the dark counts curve, or between 45 to 50 mV. 

Figures 18, 19, and 20 show the histograms of the thermal noise counts, the 
total dark counts, and the total counts under 1.0 pW received optical power 
when the threshold level was set to 47.1 mv. The counting interval was set to 1 
ms and the sample size was 50,000. The data shown on Figures 18, 19, and 20 
are subject to the quantization error of the counter used which only output 
integers. Theoretically, the numbers of counts should all follow Poisson distri- 
butions and the ratio of the variance to the mean should be equal to unity. It 
was difficult to tell from Figure 18 if the circuit thermal noise counts followed a 
Poisson distribution since there were too few counts in the 1 ms counting inter- 
val. The mean and the standard deviation measured over 1 second were about 
900 and 100, respectively. Therefore the distribution of the circuit thermal noise 
counts was clearly not Poisson. This may be attributed to the fact that the cir- 
cuit noise was somewhat correlated and the power spectrum was not perfectly 
flat, as shown in Figure 14. The distributions of the dark counts and total 
counts under 1.0 pW received optical power appeared to be Poisson within the 
quantization error of the counter, as shown in Figures 19 and 20. 

5.4. Interarrival Times of the Photon Counts 

The interarrival times of the observed counts were measured using the 
apparatus shown in Figure 21. The comparator output was divided into both 
Channels A and B of the counter which measured the time interval started by a 
pulse at Channel A and stopped by the next pulse at Channel B. The input at 
Channel A was delayed by 4.66 ns through a coax cable to avoid the counter 
being started and stopped by the same pulse. The two channels of the counter 
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had to be well balanced so that the pulse sequence seen by Channel A was the 
same as that seen by Channel B. 

Figure 22 shows a histogram of the measured interarrival time of the 
observed total number of counts under 1.0 pW received optical power. The sam- 
ple size in this measurement was 2X10 5 . The heavy solid curve in Figure 22 
represents the theoretical values of the exponential distribution of the same 
mean. The measured interarrival times fit very closely to an exponential distri- 
bution. This shows that the performance of the photon counter being tested was 
very close to that of an ideal photon counter and the detected number of counts 
follows a Poisson distribution. Figure 23 is a magnification of Figure 22 from 0 
to 100 ns. The sharp decrease near the origin was due to the dead time of the 
photon counter. The fluctuations in the histogram of Figure 22 were due to the 
finite sample size of the measurements. The average dead time can be approxi- 
mated as 15 ns after considering the coax cable delay (4.66ns) at Channel A of 
the counter. Since the dead time was a random number as mentioned in the 
previous section, the measured interarrival times could occasionally be shorter 
than the average dead time, as shown in Figure 23. 

5.5. Dynamic Range of the Photon Counter 

Table 1 shows the mean and standard deviation of the total number of 
observed counts as a function of the received optical power. The fourth column 
shows the ratio of the variance to the mean of the measured counts (excluding 
the effect of the thermal noise counts). It is seen that the detected counts had a 
near Poisson behavior up to 30 pW received optical power. The circuit seemed to 
start to saturate as the input optical power exceeded 30 pW. The threshold vol- 
tage was set to 47.1 mV so that the total number of dark counts was about 
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6~7/ms. Table 1 also lists the raw and the dead time corrected equivalent 
quantum efficiencies (excluding the dark counts). It is shown that the photon 
counter had a nearly constant equivalent quantum efficiency of about 5.0% for 
up to 10 pW received optical power or a count rate of about 2000 cts/ms 
(2,000,000 cts/sec). The total number of counts versus the received optical 
power, the data in the first and second column of Table 1, are also plotted in 
Figure 24. The curve tends to level as the received optical power becomes 
extremely low because dark counts start to dominate. The curve also starts to 
level off as the received optical power becomes relatively high because of the 
saturation of the circuit. The center segment of the curve is approximately 
linear and corresponds to an equivalent quantum efficiency of about 5.0%. This 
closely agreed with the theoretically predicted value, 6.0% ( P disc = 7.8%), accord- 
ing to (20) — (22) assuming that the average APD gain was <7=1000, APD quan- 
tum efficiency, rj = 77%, and the effective threshold current level, 

I tkre = 4.5Xcr /amp . The theoretical model given in Section 4 assumes that the fre- 
quency response of the entire circuit up to the comparator is the same as that of 
a RC lowpass filter but the experimental system contained a higher order 
lowpass filter which had a sharp cutoff frequency. Therefore the actual ampli- 
tude of photocurrent pulses might have been smaller than that given by (19). 
This is believed to be the reason why the theoretically predicted effective quan- 
tum efficiency using (19)-(22) is slightly lower than that of the experimental sys- 


tem. 
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6. Conclusions 

We have demonstrated that an APD can be used to count photons while 
being biased below its break-down point. The major advantages of this type of 
APD photon counter versus a conventional Geiger mode APD biased above the 
break-down point are the significantly shorter dead times and the total elimina- 
tion of after pulsing. A theoretical model has been developed to fully character- 
ize this type of nonbreak-down APD photon counter. The theoretically 
predicted performance agreed with the experimental measurements. The experi- 
mental nonbreak-down APD photon counter achieved an effective quantum 
efficiency of 5.0% at X = 820 nm under a dark count rate of about 7000/second. 
The counting statistics appeared to follow a Poisson distribution and the 
interarrival times fit very tightly to an exponential distribution. The measured 
dead time of the photon counter was about 15 ns and there was no after pulsing 
observed. The experimental nonbreak-down APD photon counter has already 
outperformed a typical conventional break-down type of APD photon counter. 

The performance of the experimental nonbreak-down photon counter may 
be further improved just by using more advanced commercially available elec- 
tronics components. One example is to use a higher speed comparator such as a 
GaAs device so that the full bandwidth of the APD and preamplifier (700 MHz) 
can be used. It is assumed that the equivalent threshold current is still 4.5 times 
the rms noise current so that the dark counts rate is about 7000/second. The 
APD preamplifier has to be well packaged to eliminate any EMI noise so that 
the equivalent input noise current meets the specification of the preamplifier 
(6 pA /y/Hz). The predicted effective quantum efficiency of the photon counter 
under these condition will be ^=9.2% at X=820 nm (rj Q = 77%) and 
rjQ f; = 0.72% at X= 1.06 //m (rj= 6%) with a dead time of a few nanoseconds. 
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The linear dynamic range of the photon counter should also be considerably 
larger due to the shorter dead time. Such a photon counter will be comparable 
to PMTs at X = 800 nm and outperform PMTs at X=1060nm. 

The transimpedance APD preamplifier may also be replaced with one which 
has a higher feedback resistance, for example, GigaBit Logic 16G072 
(R f =S0KQ,B 3dB =600MHz) [24]. The predicted value of P disc under these 
condition will be 24.3% which corresponds to an effective photon counting quan- 
tum efficiency of VQ t/l = 18.7% at X=820 nm and 1.46% at X= 1.06 ^/m 

under the same dark count rate, assuming the average equivalent noise tempera- 
ture of the preamplifier is T e = 1500 'I(. This photon counter will outperform 
any PMT at X > 800 nm and any APD in analog mode with excess noise factor 
F >4.1. 

If we use a state-of-the-art silicon APD [25] in which = 0.005, G = 
1000, and t] Q = 90% at X=820, the predicted value of the discrimination proba- 
bility will be P d i sc — 30.7%. The corresponding effective photon counting quan- 
tum efficiency will be VQ, lf = 27.6% at X=820nm and VQ t{f = 2.15% at 
X=1060nm at a dark count rate of less than 10, 000 /sec. This photon counter, if 
built, will clearly be superior to any PMT and APD in analog mode at the 
wavelength range. 
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Figuire 3. Avalanche break-down APD photon counter. 
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Figure 4. Normalized bias of the ML estimator, E {n }, as a function of the aver- 
age number of absorbed photons, n . 
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Figure 5. Normalized SNR of the ML estimator SNR{h ) as a function of the 
average number of absorbed photons, n. for a given dead time of r ^ and a count- 
ing interval of T seconds. 
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Figure 6. Circuit diagram of nonbreak-down APD photon counter. 




Figure 7. Photocurrent pulse shape. The dashed 
line represents the ideal exponential pulseshape 
and the solid line represents a more realistic 
pulseshape. 
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Figure 8. Probability of APD gain. 




Threshold 


Figuire 10. Average numbers of detected photons and 
dark counts as a function of the threshold. 
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Figure 11. 


Optical Setup of the Experiment. 
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Figure 12. Non break-down APD photon counter circuit 
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Figure 13. APD/Preamp DC supply. 
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Figure 14 . Power spectrum of the total circuit thermal noise measured at the 
input of the comparator. 
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Figure 16. Average waveforms of input pulses (negative going) and output pulses 
(positive going) of the comparator. 
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Figure 17. Average number and standard deviation (solid curve) of the detected 
counts in 1 ms under 1 pW received optical power as a function of the threshold 
voltage applied at the comparator. The dot-dashed and the dotted curves 
represent the circuit thermal noise counts (APG gain G=0) and the total dark 
counts (G~1000), respectively. 





Figure 18. Histogram of circuit thermal noise noise counts over a 1 ms counting 
interval. The sample size was 5X10 4 . 
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Figure 19. Histogram of the number of the total measured dark counts over a 1 
ms counting interval. The sample size was 5X10 4 . 




Figure 20. Histogram of the number of the detected counts under 1 pW received 
optical power over a 1 ms counting interval. The sample size was 5X10 4 . 




Figure 21. Apparatus used for measuring interarrival times. 
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Figure 22. Histogram of the interarrival times under 1 pW received optical power. 
The solid curve corresponds to the exponential distribution with the same mean 
(i.e. 3308 exp [— t /4.78X10 -6 ]). The sample size was 2X1 0 3 . 




Table 1. Measured counts versus received optical power. 


P(pW) 

cts/ms 

a 

var/mean* 

yj** *** 

-VI * * * 

tcorr 

0.0 

6 

3 

1.00 

N/A 

N/A 

0.1 

28 

6 

1.18 

5.1% 

5.1% 

0.3 

69 

8 

0.88 

5.0% 

5.0% 

1.0 

216 

15 

1.02 

5.2% 

5.2% 

3.0 

629 

26 

1.07 

5.1% 

5.1% 

10 

1990 

43 

0.93 

4.8% 

5.0% 

30 

5270 

72 

0.98 

4.2% 

4.6% 

40 

6510 

87 

1.16 

3.9% 

4.3% 

50 

7500 

99 

1.31 

3.6% 

4.1% 

60 

8670 

107 

1.32 

3.5% 

4.0% 

70 

9630 

127 

1.67 

3.3% 

3.9% 

100 

12.3X10 3 

153 

1.91 

3.0% 

3.6% 

200 

17.7X10 3 

214 

2.58 

2.1% 

2.9% 

300 

20.7X10 3 

261 

3.28 

1.7% 

2.4% 

400 

22.7X10 3 

294 

3.80 

1.4% 

2.1% 

500 

24.7X1 0 3 

232 

2.18 

1.2% 

1.9% 


* var/mean = {<7^—0tiur)/i mean ~ mean tkcr) with the mean and the standard deviation of 
the circuit thermal noise counts, mean t i ier =l and a th „=2, respectively. 

** r}=(mcan-mean dari )/(P in T/hf ) with the counting interval, T=lms, the received opti- 
cal power, P in = 1.0 pW, and the photon energy, A/=1.51eV. 

*** Va>TT'cud=n/[l-{ T d/T)n c ] with r d the dead time (15ns) and n c the average number of 
observed counts. 
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Figure 24. Average number of detected counts of the experimental nonbreak- 
down APD photon counter versus received optical power in pW. 



